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ABSTRACT 

Seismic  noise  generated  by  surface  winds  severely  limits 
the  usefulness  of  long-period  seismograms.  A  linear  depen¬ 
dence  is  observed  between  the  recordings  of  a  long-period 
seismograph  and  a  co-located  microbarograph  in  the  period 
range  20  to  100  seconds.  It  is  possible  to  design  a  filter, 
based  on  the  least-mean-square  method,  such  that  the  effects 
of  the  wind  generated  noise  on  the  seismograph  may  be  pre¬ 
dicted  from  the  atmospheric  pressure  changes  recorded  by  a 
microbarograph  located  at  the  same  site.  The  filter  has 
been  used  to  reduce  the  noise  level  for  seismic  data  and 
hence  increase  the  S/N  ratio.  However,  only  limited  success 
has  been  achieved  and  the  efficiency  of  the  prediction  tends 
to  depend  on  the  wind  direction,  wind  speed  and  the  variance 
of  the  wind  speed.  The  turbulent  cells  associated  with  at¬ 
mospheric  pressure  changes  are  believed  to  be  broken  down 
into  smaller  cells  as  the  speed  and  variability  of  the  wind 
increase.  The  increase  of  high  frequency  energy  contributed 
by  smaller  cell  disturbances  also  increases  the  nonlinear 
part  of  the  transfer  function  and  hence  breaks  down  the  pre¬ 
diction  ability  of  the  filter. 
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The  background  noise  in  the  period  range  20  to  100 
seconds  has  always  caused  difficulties  in  the  analysis  of 
the  surface  wave  data  on  long-period  seismograms.  A  corre¬ 
lation  has  often  been  observed  between  changes  in  the  magni¬ 
tude  of  the  long-period  seismic  noise  and  the  atmospheric 
pressure  fields.  By  studying  the  disturbances  on  the  long- 
period  vertical  seismograph  and  the  atmospheric  pressure 
fluctuations  on  a  nearby  microbarograph,  Crary  and  Ewing 
(1952)  concluded  that  the  correlation  could  be  explained  by 
atmospheric  buoyancy  effects  on  the  seismometer.  Great  im¬ 
provement  in  the  performance  of  long-period  vertical  seis¬ 
mographs  were  obtained  by  elimination  of  buoyancy  effects 
(Ewing  &  Press,  1953).  Steps  have  been  taken  to  avoid  direct 
pressure  effects,  such  as  the  seismometer  being  built  in  a 
tightly  sealed,  rigid  steel  case,  which  is  then  enclosed  in 
a  sealed  steel  vault.  Still  substantial  correlation  persists 
(Haubrich  &  Mackenzie  1965,  Capon  1969)  which  is  due  to  the 
elastic  deformation  of  the  ground  caused  by  the  atmospheric 
pressure  changes. 

Noise  recorded  by  long-period  seismometers  consists 
of  both  propagating  and  non-propagating  components  (Ziolkowski, 
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1973).  Rayleigh  waves  generated  by  sources  such  as  the  action 
of  surf  on  the  coast  and  from  cultural  activity  are  propa¬ 
gating  noises.  Approximately  half  the  noise  observed  by 
long-period  seismometers  at  LAS a*  (Capon,  1969)  is  non-propa¬ 
gating;  that  is,  it  is  incoherent  over  distances  of  a  few 
kilometers  and  is  not  of  cultural  or  oceanic  origin. 

Previous  studies  have  included  mathematical  models 
which  describe  the  non-propagating  component  of  the  noise 
field.  A  plane  pressure  wave  is  assumed  (Sorrells,  1971)  to 
propagate  at  the  mean  wind  speed  over  different  models  of 
layered  half-spaces.  These  numerical  studies  indicate  that 
pressure  waves  can  contribute  significantly  to  the  output  of 
the  long-period  vertical  seismograph,  and  also  create  severe 
tilt  noise  on  long-period  horizontal  seismographs.  The 
pressure  response  of  the  earth  tends  to  increase  monotonically 
as  the  wavelength  of  the  pressure  oscillation  increases.  The 
seismic  disturbances  created  by  plane  pressure  waves  decay 
rapidly  with  depth.  Long-period  seismographs  have  been 
situated  both  on  the  surface  of  the  earth  and  in  a  mine  at  a 
depth  of  183  meters  at  Grand  Saline,  Texas,  by  Sorrells  et  ad 
(1971).  The  observations  showed  that  during  windy  intervals, 
the  local  pressure  change  is  strongly  correlated  with  the 

*  Large  Aperture  Seismic  Array  (LASA) ,  Billings,  Montana 
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noise  recorded  on  the  surface  and  is  incoherent  with  the  noise 
recorded  in  the  mine  for  the  period  range  20  to  100  seconds 
(Fig.  1).  The  ground  noise  on  the  long-period  seismogram 
is  thus  of  local  origin,  non-propagating  and  largely  confined 
to  the  surface  of  the  earth.  Similar  results  have  been  ob¬ 
tained  by  Savino  et.  al.  (1972)  at  a  mine  near  Ogdensburg,  New 
Jersey.  The  noise  recorded  in  the  mine  remained  a  continuous 
low-level  disturbance  while  the  noise  level  on  the  surface 
varied  by  about  10  dB  over  a  few  days  depending  on  the  meteor¬ 
ological  conditions.  The  wind  generated  noise  attenuates 
rapidly  with  depth  and  can  be  virtually  eliminated  by  install¬ 
ing  the  seismographs  at  a  depth  of  several  hundred  meters 
(Sorrells  1971). 

The  coherence  calculated  from  the  recordings  of  a  long- 
period  seismograph  and  a  nearby  microbarograph  indicates  that 
there  is  a  linearly  dependent  component  between  the  two  time 
series;  therefore,  it  is  possible  to  estimate  the  pressure 
generated  earth  noise  by  applying  a  linear  transfer  function 
to  the  microbarograph  output.  The  transfer  function,  called 
a  Wiener  filter,  is  based  on  the  least-mean-square  method 
(Burg  1964) .  It  has  been  applied  to  long-period  seismic 
data  by  Ziolkowski  (1973)  and  Douze  ei  aL.  (1975).  Indiscrim¬ 
inant  applications  of  the  Wiener  filter  by  these  authors  re¬ 
sulted  in  inconsistent  results.  In  some  cases,  dramatic 


P*' *00  (a) 


Fig.  1  Data  recorded  during  a  windy  period  with 

mean  wind  speed  7.8  m/s. 

(a)  Power  spectrum  of  the  noise  recorded  by  base 
microbarograph . 

(b)  A  comparison  of  the  power  spectral  density 
estimates  of  the  noise  recorded  by  surface  and 
mine  vertical  seismographs. 

(c)  Low  coherence  between  surface  and  mine  vertical 
seismographs. 

(d)  Base' microbarograph  is  highly  correlated  to 
surface  vertical  seismograph  and  there  is  no 
coherence  between  base  microbarograph  and  mine 
record.  (from  Sorrells  ejt  al.  1971) 
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improvements  in  signal-to-noise  ratio  was  achieved?  in  other 
cases,  no  improvements  whatsoever  resulted.  It  is  the  pur¬ 
pose  of  this  study  to  identify  the  cause  or  causes  of  the 
inconsistent  results,  and  to  determine  the  optimum  manner 
with  which  Wiener  prediction  filtering  can  be  used  with 
state-of-the-art  long-period  seismic  recording. 


NOISE  PREDICTION  TECHNIQUES 


The  noise  caused  by  wind  generated  pressure  changes 
recorded  on  long-period  seismographs  may  be  highly  coherent 
with  the  microbarograph  recordings  at  the  same  site  when 
the  wind  blows  strongly.  The  linear  dependence  between  the 
outputs  of  the  two  instruments  leads  to  the  development  of 
a  filter  which  predicts  the  wind  related  noise  in  the  seismic 
data.  A  Wiener  filter,  also  known  as  an  optimum  shaping 
filter,  is  designed  to  produce  a  certain  desired  output  from 
some  input  signal.  In  the  case  of  this  experiment,  the  noise 
data  on  the  seismogram  is  the  desired  output  for  a  shaping 
filter  while  the  microbarograph  data  is  the  input  signal. 

The  filter  is  a  linear  system  which  best  transfers  the  micro¬ 
barograph  signal  to  ground  noise  under  the  least  mean  square 
criterion.  The  mathematical  operation  is  described  in  the 
following  sections . 

The  output  of  the  seismograph  S (t)  can  be  divided  into 
two  parts,  N{t)  and  E(t).  N{t)  is  the  time  series  which  is 
linearly  related  to  the  microbarograph  output  P(t),  and  is 
caused  by  local  pressure  fluctuations  in  the  atmosphere.  On 
the  other  hand  E(t)  is  not  linearly  related  to  the  pressure 
and  contains  noises  generated  from  other  sources  and  the 


desired  signals,  thus 

S  (t )  =  N(t)  +'  E  ( t )  (Eq.  1) 

Since  N(t)  is  experimentally  correlated  with  P(t),  a  Wiener 
filter  F(t)  can  be  determined  by  taking  N(t)  to  be  the  de¬ 
sired  output  with  input  P(t);  thus  P(t)  convolved  with  the 
filter  yields  N(t),  an  estimate  of  N(t). 

P(t)  *  F (t)  =  N(t)  (Eq.  2) 1 

To  design  the  filter  F(t),  a  time  period  should  be 
chosen  from  intervals  with  high  wind  speed  and  no  other 
obvious  seismic  signals  so  that  most  of  the  noise  energy 
on  the  long-period  seismogram  is  wind  related;  i.e.  E(t)  = 
The  seismic  output  S (t)  in  this  time  period  is  then  approxi 
mately  equal  to  wind  noise  N(t),  the  desired  output  for  the 
shaping  filter.  Equations  1  and  2  can  be  written  as 
S(t)  ~N(t)  and  P(t)  *  F (t )  =  S(t) 

Steps  to  generate  the  filter  are  described  in  detail  by 
Kanasewich  (1973)  and  Robinson  (1967),  and  are  summarized 


as  follows. 

Let  microbarograph  P(t)  =  Input  t  =  l, . ,T 

seismograph  S(t)  =  desired  output  t  =  1, . ,T 

f(k)  =  unknown  filter  coefficients  k  =  1, . ,K 

A(t)  *  output  after  filtering  t  =  1, . ,T 


1.  The  symbol  *  indicates  convolution. 
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The  actual  output  A(t)  is  the  convolution  of  microbarograph 
P(t)  and  the  unknown  filter  f(k),  that  is 
P(t)  *  f(t)  =  A (t ) 

K 

A ( t )  =  Z  f(k)  P(t-k)  t  =  1, . . .T+K-l 

k=l 

The  error  between  A(t)  and  desired  output  S(t)  is  defined 
as 

e(t)  =  A (t )  -  S(t) 

e(t)  =  £  Z  f(k)  P(t-k)J  -  S(t)  t*l _ T+K-l 

The  errors  are  squared  and  summed  for  the  time  series  as  E. 


e(t)2  = 

f  E 

f(k)  p(t-k)  1 

-S(t) 

Lk=l 

J 

T+K-l 

E  =  Z  e<t) 

t=L 


T+K-l  ( p  K 

E  -  Z  Z  fOO  P(t 
t=l  I  L  k=l 

The  filter  is  best  fitted  with  minimum  error  if  the  partial 
derivatives  of  E.  with  respect  to  each  of  the  filter  coeffi¬ 
cients  are  equated  to  zero. 
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dE 


df(j) 


T+K-1 

rK  i 

=  2  2 

2  f(k)  P(t-k) 

-S(t)  P(t-j)  =  0 

t=l 

L  k=l  J 

which  gives  K 

2  f(k) 
k=l 


r  t+k-1 

2 

.  t=l 


P(t-k)  P(t 


]  T+K-1 

-  2 

t=l 


S(t)  P(t-j) 


j  -  1 . K  (Eq.  3) 


T+K-1 

where  2  P(t-k)  P(t-j)  =  <t>  is  the  autocorrelation  of 
t=l  j-* 

P(t)  and 


T+K-1 

2  S (t)  P(t-j)  =  C_j  is  the  crosscorrelation  of 
t=l 

S (t)  &  P(t) . 

By  substituting  4>  and  C_j  into  Eq.  3,  we  have 
K 

2  f(k)  4>(j-k)  =  C(-j)  j=l, - -  K  (Eq.  4). 

k=l 


Hence  the  desired  filter  coefficients  are  calculated  by 
using  the  autocorrelation  of  microbarograph  and  the  cross¬ 
correlation  of  microbarograph  and  seismic  noise  data.  Eq.  4 
represents  a  set  of  K  simultaneous  equations  and  can  be 
written  out  in  the  matrix  form, 


’0(0)  0(1)  ... 

♦d)  0(0)  ... 

. .  0 (k-1) * 

•  • 

¥ 

0(k-l)  _ 

0(0) 

m 

A. 

.C*_ 

The  K  filter  coefficients  are  unknown  and  can  be  calculated 

by  inverting  the  KxK  matrix,  but  the  computation  work  is  great 

for  a  large  number  of  filter  coefficients.  The  problem  is 

simplified  by  using  the  recursive  method  first  introduced 

by  Levinson  (1947),  in  which  one  solves  the  equations  in  a 

step-wise  manner.  We  start  by  computing  a  one-point  filter, 

and  then  the  steps  for  filters  with  2,  3,  _ _  K  points  are 

calculated  successively.  The  machine  time  required  to  solve 

the  simultaneous  equations  for  a  filter  with  K  coefficients 

2 

is  proportional  to  K  using  the  recursive  method,  as  compared 
to  using  the  conventional  method.  The  computer  storage 
space  is  also  condensed  by  using  this  recursive  method. 

An  infinitely  long  filter  is  required  to  produce  exactly 
a  desired  output  and  it  follows  that  E  decreases  as  filter 
length  increases.  The  error  of  the  filter  is  computed  at 
each  filter  length  during  the  iteration,  and  usually  the 
error  will  decrease  rapidly  and  then  level,  off.  The  least 
number  of  filter  points  required  for  minimum  error  can  be 
determined  by  experiment. 

The  Wiener  filter  is  then  applied  to  a  time  period  of 
interest  which  is  within  two  or  three  hours  of  the  time  period 
from  which  the  filter  was  obtained.  The  wind-related  noise 
N(t)  can  be  estimated  by  convolving  the  microbarograph  output 
with  the  filter  to  obtain  N(t).  Subtracting  N(t)  from  the 


seismic  data  S(t)  yields  an  improvement  in  signal -to-noise 
ratio,  providing  that  the  linearly  dependent  portion  of  the 
noise  field  is  significant. 


EXPERIMENTS  AND  RESULTS 

The  site  of  the  experiments  was  McKinney,  Texas, 
about  30  miles  north  of  Dallas.  At  McKinney,  seismic  and 
meteorological  data  were  recorded  continuously  on  thirty-two 
channels  of  digital,  magnetic  tape.  Channels  containing 
data  from  the  long-period  vertical,  long-period  E-W,  long- 
period  N-S  seismographs;  and  the  microbarograph  and  wind 
speed  and  wind  direction  indicators  were  used  in  this  experi¬ 
ment.  The  wind  speed  in  the  McKinney  area  during  the  year 
1974  was  generally  low.  Prolonged  periods  of  wind  speeds  in 
excess  of  5  meters/second  (about  11  miles/hr)  were  rare. 

Most  of  the  earthquakes  recorded  at  McKinney  with 
magnitudes  (m*,)  greater  than  4.5  show  good  Rayleigh  wave 
dispersion  patterns  on  the  long-period  vertical  component 
(Fig.  2),  but  well  dispersed  Love  waves  are  much  more  diffi¬ 
cult  to  observe.  Surface  horizontal  seismographs  are  much 
more  susceptible  to  wind-generated  noise  than  are  vertical 
instruments  and  as  a  result  generally  exhibit  a  much  higher 
noise  level.  The  flat  dispersion  curve  within  the  pass 
band  of  the  instrument  for  oceanic  Love  waves  results  in  a 
poorly  dispersed  wave  which  does  not  contrast  greatly  with 
the  noise.  Nevertheless,  several  examples  were  obtained. 
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As  noted  previously,  the  noise  related  to  the  atmospheric 
effects  is  dominated  by  quasi-static  deformations  in  response 
to  wind  generated  pressure  changes  during  windy  intervals 
with  moderately  high  wind  speed  (>4  m/s)  (Douze  e£.  aJL,  1975). 
The  wind  moving  over  the  surface  also  causes  tilts  on  the 
horizontal  components  which  depend  on  the  wind  direction. 

The  Preliminary  Determination  of  Epicenters  (PDE)  of  the 
USGS  and  the  meteorological  data  for  McKinney  during  the  year 
1974  were  reviewed.  Seventeen  earthquakes  were  selected 
which  were  recorded  during  different  conditions  of  wind  speed 
and  variability.  Suitable  filters  were  designed  and  applied 
to  each  event.  The  details  of  the  filter  design  and  the 
data  processing  is  described  in  Fig.  3.  The  magnitudes  (m^) 
of  the  events  selected  ranged  between  5.0  and  6.0.  A  period 
of  noise  data  within  three  hours  before  or  after  the  signal 
with  similar  wind  conditions  was  selected  for  filter  design. 
PDE  sheets  were  checked  to  assure  that  no  small  events  were 
included  in  the  noise  sample.  The  3-component  data  were 
rotated  for  each  event  so  that  the  horizontal  seismographs 
were  effectively  oriented  along  directions  radial  and  trans¬ 
verse  to  the  travel  path.  Both  earthquake  and  noise  data 
were  rotated  and  a  Wiener  filter  was  designed  for  each  channel 
using  the  microbarograph  data  as  an  input.  The  output  of  the 
filter  design  program  shows  the  error  of  the  prediction.  An 
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error  less  than  0.50,  which  means  more  than  50%  of  the  noise 
energy  is  predicted,  is  considered  to  be  acceptable.  Table 
1  shows  some  of  the  results  for  earthquakes  using  Wiener 
filtering.  The  errors  for  all  three  rotated  components  are 
included.  Event  0061  which  is  an  earthquake  in  the  Solomon 
Islands  with  magnitude  5.2  has  an  acceptably  low  prediction 
error,  0.178  on  the  transverse  direction.  Fig.  4  shows  the  seis¬ 
mograms  before  and  after  Wiener  filtering.  There  is  a  large 
reduction  in  the  noise  level,  and  the  presence  of  a  Love  wave 
signal  can  be  recognized  on  the  filtered  transverse  record. 

All  the  other  examples  show  poorer  prediction  of  the 
noise  energy  on  one  or  more  channels.  Several  reasons  for 
the  failure  of  this  approach  are  suggested. 

1.  Douze  and  Sorrells  (1975)  found  that  the  prediction 
capability  increased  with  the  mean  wind  speed  and  the  Wiener 
filter  was  most  effective  when  the  wind  speed  was  greater 
than  5  m/s . 

2.  When  high  wind  velocity  prevails,  the  variability 
of  the  wind  speed  and  direction  may  reduce  filter  effective¬ 
ness. 

An  effort  was  made  to  determine  if  the  above  suggestions 


were  valid,  and  if  so,  to  determine  the  conditions  of  wind 
•peed  and  variability  where  Weiner  filtering  can  be  expected 
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Table  1.  Prediction  errors  on  three  rotated  channels 
for  earthquake  data 


Event 

Day-time 

Angle  (Deg) 
rotated 

Wind 

dir.  speed 

Prediction 
V  T 

error 

R 

0008 

118-1930 

275 

S 

(M/S) 

6.5 

.7  54 

.740 

.898 

0015 

116-2250 

11.3 

S 

5.8 

.739 

.789 

.649 

0019 

119-1417 

301 

S 

6.9 

.523 

.615 

.691 

0050 

131-1653 

302 

NE 

4.2 

.767 

.648 

.858 

0061 

136-0630 

268 

S 

6.2 

.352 

.178 

.487 

0068 

255-1100 

155 

S 

4.5 

.678 

.325 

.228 

0069 

264-0900 

277 

NE 

4.5 

.413 

.251 

.314 

0070 

264-2200 

150 

NE 

4.5 

.336 

.465 

.582 

0071 

265-1500 

304 

E 

4.5 

.337 

.831 

.809 

0072 

151-0300 

2.9 

S 

5.6 

.341 

.487 

.241 

0073 

118-1240 

21.4 

S 

6.0 

.552 

.763 

.626 

0074 

118-1640 

251 

s 

5.6 

.848 

.892 

.941 

0075 

118-2300 

305 

s 

6.8 

.682 

.870 

.856 

0076 

151-0210 

314 

s 

5.6 

.429 

.389 

.318 
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to  be  effective.  Additional  noise  data  were  collected  and 
the  filter  along  with  the  prediction  error  for  each  channel 
was  generated.  Again,  each  noise  sample  was  carefully  chosen 
to  exclude  any  identifiable  small  signals,  propagating  noise, 
or  instrumental  malfunction.  Table  2  shows  the  result  of 
the  experiment,  including  the  wind  speed  and  its  standard  de¬ 
viation.  All  the  noise  data  used  for  prediction  are  1024 
seconds  long,  and  the  length  of  the  filter  is  232  points. 

Among  most  of  the  examples,  the  prediction  on  the  horizontal 
components  is  better  than  that  of  the  vertical  component. 

This  was  not  unexpected  because  horizontal  seismographs  are 
more  susceptible  to  wind  related  earth  tilts  (Sorrells,  1971; 
Sorrells  and  Goforth,  1973)  which  cause  larger  apparent  dis¬ 
placements  at  periods  greater  than  20  seconds  than  are  seen 
on  the  vertical  component.  The  tilting  effects  caused  by 
the  pressure  changes  also  depend  on  the  wind  direction  (Sorrells 
and  Goforth,  1973).  It  can  be  seen  from  Table  2  that  lower 
prediction  errors  for  the  horizontal  channels  always  occur 
with  the  instruments  aligned  more  nearly  in  the  direction  of 
the  wind.  An  example  of  noise  reduction  is  shown  in  Fig.  5. 

Fig.  6  shows  the  plot  of  the  percentage  of  total  energy 
of  the  noise  that  can  be  predicted  on  the  horizontal  channels 
versus  the  average  wind  velocity.  The  percentage  of  noise 
prediction  is  calculated  by  one  minus  the  average  prediction 

\ 
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Table  2.  Wind  speeds,  wind  speed  deviations  and 
prediction  errors  for  noise  data 


Event 

Start 

Angle  (Deg) 
Rotated 

Wind 

Soeed 

(M/S) 

Prediction 

Errors 

No. 

Day-time 

Dir. 

Ave . 

Dev . 

Ch . 9  Chll 

Z  H 

Chi  2 

H 

01 

121-0900 

148 

S 

0.47 

0.35 

.875  ND 

.875 

02 

125-0417 

258 

NE 

3.25 

0.57 

.701  .438 

.489 

03 

126-0738' 

250 

N 

1.00 

0.22 

.825  .760 

.779 

04 

127-0954 

149 

N 

0.41 

0.50 

.794  .692 

.559 

05 

131-1653 

302 

NE 

2.92 

0.71 

.767  .648 

.858 

06 

136-0630 

268 

S 

6.94 

1.20 

.352  .178 

.487 

07 

140-2233 

310 

S 

5.46 

0.94 

.558  .566 

.396 

08 

151-0210 

314 

s 

5.88 

0.83 

.429  .389 

.318 

90 

151-0300 

2.9 

s 

5.97 

0.74 

.341  .487 

.241 

10 

255-1100 

155 

s 

3.27 

0.37 

.678  .325 

.228 

11 

264-0900 

277 

NE 

3.78 

0.68 

.413  .251 

.314 

12 

264-2200 

150 

NE 

3.68 

0.98 

.336  .465 

.582 

13 

51-0700 

S 

6.82 

0.99 

.286  .221 

.166 

14 

51-1300 

S 

6.58 

0.99 

.457  .279 

.148 

15 

77-0500 

s 

5.77 

0.60 

.280  .251 

.178 

16 

77-0600 

s 

5.76 

0.57 

.323  .241 

.293 

17 

136-0425 

s 

7.03 

0.84 

.573  .608? 

.175 

* 
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of  all  three  figures  are  in  the  same  scale.  The 
wind  was  6.58  m/s  from  the  south.  The  prediction 
error  on  N-S  direction  is  n.148. 


PERCENTAGE  OF  NOISE  PREDICTED 
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error  for  the  two  horizontal  instruments.  The  predictions  are 
scattered  and  the  scatter  increases  with  increasing  wind  speed. 
The  scattering* -effect  is  shown  to  be  caused  by  the  variability 
of  the  wind  speed  about  the  mean  in  Fig.  7.  Fig.  7  is  a  plot 
of  the  percentage  of  prediction  versus  the  standard  deviation 
of  mean  wind  velocity.  The  numbers  represent  the  mean  wind 
speed  for  each  event.  The  dashed  lines  are  connected  for 
points  with  wind  speed  of  3-4  m/s,  5-6  m/s,  and  6-7  m/s  re¬ 
spectively.  It  seems  that  the  wind  noise  prediction  from 
Wiener  filtering  is  not  only  a  function  of  the  mean  wind  speed, 
but  also  a  function  of  the  deviation  about  the  mean  of  the  wind 
speed.  The  prediction  error  is  small  for  each  velocity  cate¬ 
gory  for  a  low  variance  of  velocity  but  the  error  increases 
rapidly  with  increasing  deviation  in  velocity.  The  higher  the 
wind  speed,  the  larger  the  deviation  of  wind  speed  permissable 
for  successful  Wiener  filtering. 


PERCENTAGE  OF  NOISE  PREDICTED 
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Fig.  7.  The  percentage  of  noise  energy  predicted  on  horizontal 
components  versus  the  standard  deviation  of  mean  wind 
speed.  The  numbers  in  the  figure  show  the  mean  wind 
velocity  in  meters/second  for  each  event.  Dash  lines 
are  connected  by  points  of  velocity  3-4m/s,  5-6  m/s, 
and  6-7  ro/s  respectively. 


DISCUSSION 


From  the  experimental  data  of  J.  McDonald  and  E.  Herrin 
(1975),  it  is  shown  that  there  is  a  relationship  between  the 
measured  wind  speed  and  the  associated  micropressure  level 
(Fig.  8).  The  pressure  energy  increases  as  the  wind  speed 
increases,  and  the  resulting  non-propagation  contribution  to 
the  long-period  seismic  noise  field  constitutes  a  greater 
percentage  of  the  total  noise  field.  Thus,  in  general,  the 
harder  the  wind  blows,  the  more  effective  Weiner  filtering 
will  be. 

The  wind,  which  is  the  motion  of  the  air  in  the  bound¬ 
ary  layer  of  the  earth  caused  by  the  atmospheric  pressure 
changes,  is  turbulent  and  consists  of  turbulent  cells  of  all 
dimensions  (J.A.  McDonald  et  al .  1971)  .  The  bigger  cells 
cause  seismic  background  noise  in  the  longer  period  range, 
while  the  smaller  cells  cause  noise  in  the  shorter  period 
range.  The  turbulence  has  a  wide  range  of  subsidiary  motions 
with  different  scale  lengths  and  the  energy  is  continuously 
passing  from  the  larger  to  the  smaller  scale  lengths.  A 
large  scale  motion  in  the  atmosphere  is  unstable  and  breaks 
down  to  motions  of  smaller  scale  length;  that  is,  to  smaller 
eddies . 
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Period  (s) 
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In  Fig.  9,  the  power  spectra  of  micropressures  recorded 
during  different  wind  conditions  have  been  plotted.  The  wind 
speeds  are  approximately  5.7  m/s  for  figure  (a)  and  6.5  m/s 
for  figure  (b) .  In  both  plots,  the  high  frequency  component 
of  the  energy  is  most  affected  by  a  high  variance  of  wind 
speed,  B  and  D.  The  increase  of  high  frequency  energy  reduces 
the  fraction  of  the  energy  in  the  pressure  field  which  is 
linearly  related  to  the  seismic  noise.  The  high  frequency 
disturbances  associated  with  the  small  scale  motions  in  the 
atmosphere  thus  interfere  with  the  prediction  of  ground  noise. 

The  linear  coherence  of  the  energies  of  two  time  series 
indicates  that  part  of  one  time  series  is  linearly  dependent 
on  the  other  time  series  while  the  whole  process  relating 
them  might  be  a  non-linear  one  (Ziolkowski,  1973).  In  other 
words,  the  linear  component  of  one  time  series  can  be  transformed 
from  the  other  with  a  linear  time-domain  filter.  Tie  design 
of  the  Wiener  filter  is  based  on  the  linear  dependence  between 
the  ground  noise  and  microbarograph  output  in  the  period  range 
20  to  100  seconds.  In  a  situation  with  variable  wind,  the  high 
frequency  energy  will  dominate  the  total  energy  spectrum  of 
the  microbarograph  (Fig.  9).  The  increase  of  high  frequency 
energy,  which  is  also  the  incoherent  part  of  the  transform 
process,  degrades  the  linear  relation  between  the  two  time 
series,  and  hence  causes  increasing  prediction  error  with 
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increasing  variance  o£  wind  speed. 

Now  let  us  take  a  look  at  the  mathematical  representation 
of  the  process.  For  any  two  time  series,  a  function  can  be 
found  such  that 

S(t)  =  Function  (P(t)) 

where  S(t)  is  the  seismic  noise  data  and  P(t)  is  the  micro¬ 
barograph  data  as  we  used  in  the  previous  section.  The 
function  can  be  expanded  in  a  polynomial  form, 

S(t)  =  a1P(t)  +  a2P(t)2  +  a3P(t)3  +  . 

=  a^Pft)  +  higher  order  terms 
=  fx(t)  +  f2(t) 

with  f^(t)  =  a^Ptt)  which  is  linearly  related  to  S ( t )  and 
f2(t)  are  the  higher  order  terms  which  are  not  linearly  re¬ 
lated  to  S (t) . 

After  taking  Fourier  transforms,  we  get 
S (f )  -  Fx(f)  +  F2(f). 

In  our  experiment,  F^f)  is  primarily  the  energy  in  the  period 
range  20  to  100  seconds  because  of  the  linear  dependence  ex¬ 
isting  between  S(t)  and  P(t)  in  this  range.  F2(t)  is  then 
the  higher  order  high  frequency  portion  for  which  Wiener 
filtering  is  ineffective.  During  the  intervals  with  steady 
wind,  F2(f)  is  small  enough  to  be  neglected.  The  linear 
transfer  function  between  the  time  series  S(t)  and  P(t)  is 
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essentially  the  linear  part  f^(t)  and  we  expect  very  good 
prediction  of  the  noise  energy  on  the  seismogram.  Once  the 


high  frequency  energy  increases  due  to  the  breaking  down  of 
the  turbulent  cells  to  small  eddies,  the  nonlinear  energy 
Fjff)  increases  (Fig.  9)  and  dominates  the  whole  process. 
The  linear  filter  is  then  working  on  a  non-linear  system 
and  thus  fails  to  predict  noise  accurately. 


CONCLUS ION 


Most  of  the  noise  recorded  on  the  long-period  seismo¬ 
graph  appears  to  be  generated  locally  by  pressure  disturbances 
in  the  atmosphere  moving  at  the  mean  wind  speed.  There  is  a 
component  of  the  seismic  noise  field  which  is  linearly  related 
to  the  atmospheric  pressure  field  in  the  period  range  20  to  100 
seconds.  An  optimum  Wiener  filter  is  designed  to  predict  and 
reduce  the  wind  noise  by  assuming  that  the  nonlinear  part  within 
the  period  range  is  small  and  can  be  neglected.  During  inter¬ 
vals  with  unsteady  meteorological  conditions  the  turbulent 
cells  are  believed  to  break  down  into  small  eddies  which  in¬ 
crease  the  amount  of  high  frequency  energy  and  lead  to  a  sig¬ 
nificant  non-linear  component  in  the  relation  between  the 
pressure  field  and  seismic  noise  in  the  period  range  20  to  100 
seconds.  Thus  the  linear  Wiener  filter  fails  in  noise  pre¬ 
diction  for  time  periods  in  which  the  wind  speed  is  highly 
variable.  Some  other  technique,  such  as  deep  burial  of  the 
seismometers  is  required  to  suppress  the  pressure  generated 
noise  under  these  conditions.  For  areas  characterized  by 
high  and/or  steady  wind,  Wiener  filtering  should  be  a  useful 
tool  for  decreasing  background  seismic  noise  on  long  period 
(20-100  sec)  surface  instruments. 
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APPENDIX 

Descriptions  of  the  program  names  used  in  Fig.  3. 

DISPLAY:  This  program  can  show  simultaneously  any  three 

seismograms  on  the  computer  scope. 

RAYTHEON:  This  program  reads  from  a  Raytheon  TC  200  Field 

tape  and  converts  the  data  to  TAPE  I/O  (1025) 
format . 

COORS:  The  program  computes  the  distance  and  azimuth 

between  reference  points. 

PLOT  TAPE 

I/O:  This  program  plots  up  to  13  blocks  of  a  channel 

on  a  data  tape  with  TAPE  I/O  (1025)  format. 
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